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ABSTRACT 

The cosmic infrared background (IRB) at wavelengths between 1 /im and 3 pro- 
vides a useful probe of early star-formation prior to and during reionization. To explain 
the high optical depth to electron scattering, as measured by the Wilkinson Microwave 
Anisotropy Probe (WMAP), one requires significant star-formation activity at red- 
shifts 10 and higher. In addition to massive stars, the IRB flux may be contributed 
by a population of early miniquasars. We study the relative contributions from first 
stars, supernovae and quasars to the IRB for reasonable star formation rates at high 
redshift. If miniquasars radiate efficiently at the Eddington-limit, current background 
measurements limit the fraction of mass in first stars that is converted to seed black 
holes to be roughly less than 10%. In the case of supernovae, though an individual su- 
pernova is much brighter than the progenitor star, due to the shorter lifetime of order 
few months, the fractional contribution to the IRB remains at a level of 10% and be- 
low when compared to the same contribution from stars. The bright supernovae may, 
however, be directly detectable by future large ground-based and space telescopes. 
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1 INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) has 
provided strong evidence for an optical depth for electron 
scattering of 0.17 ± 0.04 based on the large scale polarization 
pattern related to rescattering of Cosmic Microwave Back- 
ground (CMB) photons (Kogut et al. 2003). If the reioniza- 
tion process is described as instantaneous and homogeneous, 
the measured optical depth implies a reionization redshift 
of ~ 17 ± 5 in a spatially flat universe. Interestingly, the 
derived redshift for reionization is at the high end of expec- 
tations in models in which stellar populations are dominant 
reionization sources (e.g., Cen 2003; Fukugita & Kawasaki 
2003; Venkatesan, Tumlinson & ShuU 2003; Wyithe & Loeb 
2003; Sokasian et al. 2004). Thus, in some of these models, 
primordial stars, the so-called Population III stars, are as- 
sumed to be rather massive and hence have a large UV pho- 
ton emission rate, as suggested by recent theoretical studies 
(e.g. Abel, Bryan & Norman 2002; Bromm, Coppi & Larson 
2002; Bromm, Kudritzki & Loeb 2001; Schaerer 2002, 2003). 

The formation of such massive stars in the early uni- 
verse has many important cosmological implications (e.g.. 



Carr, Bond, & Arnett 1984; Oh, Cooray Kamionkowski 
2003; Yoshida, Bromm & Hernquist 2004). Because radi- 
ation below the Lyman-limit is absorbed by neutral hydro- 
gen, luminous sources at redshifts 10 and higher are gener- 
ally seen only in the near-IR band. In particular, if Pop III 
stars are formed primarily at redshifts between 10 and 30, 
they are expected to contribute to the infrared background 
(IRB) light at wavelengths between 1 and 5 /xm at present 
epoch (Bond, Carr & Hogan 1986). Recent estimates based 
on theoretical models suggest that a large fraction of the 
IRB total intensity may indeed be due to these stars (Santos 
et al. 2002; Salvaterra & Ferrara 2003; Cooray et al. 2004). 
A substantial IRB could arise from the Pop III stars not 
only due to the direct emission associated with these stars, 
but also due to indirect processes that lead to free-free and 
Lyman-Q emission from the ionized nebulae, or Hll regions, 
surrounding these stars. 

Observationally, about 20-40% of the near-IR flux (be- 
tween ~ 1 /im and 3 fim) has been resolved by point sources 
(Cambresy et al. 2001; Totani et al. 2001); while J- and 
K-band source counts go down to AB-magnitudes of ~ 28 
and 25, respectively, the cumulative surface brightness con- 
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verges by a magnitude around 23, say, in the K-band (see for 
a review of existing data by Pozzetti & Madau 2000). The 
missing IR flux could be either due to a high-z population of 
galaxies, such as the proto-galaxy population related to first 
stars (Santos ct al. 2002; Salvatcrra & Fcrrara 2003), or a 
population of low surface brightness gala^xies at nearby red- 
shifts. While there still remains uncertainties, one can poten- 
tially understand the presence of first luminous sources and 
some details related to their population, such as the number 
density, by characterizing the background light from optical 
to infrared at wavelengths of a few micron. In this wave- 
length range, final products of these massive stars are also 
expected to contribute to IRB, in addition to the direct emis- 
sion from stars; Radiation from supcrnovac and primordial 
black holes, in the form of miniquasars, may be significant 
sources of emission. The role of final products in contribut- 
ing to the IRB, relative to the progenitor stars, depends on 
the initial mass function of the first stars. 

As the end state of a massive star above ~ 26OM0 is a 
complete collapse to a black hole (Heger & Woosley 2002), 
if the first stellar population is extremely top-heavy with all 
mass above this limit, then one expects roughly a similar 
number density of black holes as that of stars. These black 
holes grow via accretion and/or mergers and may partly be 
responsible for early reionization if they radiate as mini- 
quasars (Haiman & Loeb 1998; Madau ct al. 2004; Ricotti 
& Ostriker 2003). The UV output from such quasars, if not 
absorbed heavily by the surrounding torii, can contribute 
both to reionization and to the IR flux as viewed today. 
The unresolved soft X-ray background at energies between 
0.5 and 2 keV may rule out a largo density of miniquasars 
at high redshifts as required to reionize the universe from 
the quasar emission alone (Dijkstra, Haiman & Loeb 2004) 
unless reionization by miniquasars, within uncertainties of 
models used in Dijkstra et al., is completed by a redshift 
of 20. Here, we suggest that if all massive stars are eventu- 
ally converted to black holes and accrete efficiently at the 
Eddington-rate (as would be the case if miniquasars are re- 
sponsible for eventual reionization of the Universe) , then the 
IR background is overproduced. A more realistic situation 
could be that the mass function is broad, and only a small 
fraction of the total initial stellar mass is converted to black 
holes; current understanding of the IRB may limit the mass 
fraction of stars above 260 Mq to a level less than 10%. 
This constraint can bo relaxed if miniquasars radiate less 
efficiently below the Eddington- limit; such a scenario may 
be viable if stellar sources dominate the UV photon produc- 
tion and are, thus, responsible for reionization, though some 
fraction of stellar mass is converted to miniquasars. 

The initial mass function of the first stars could po- 
tentially be dominated by stars with masses in the range 
between 140 and 260 M© . In this case, the final product will 
be a complete disruption of the star via the pair-instability 
process (Heger & Woosley 2002). In the case of such mas- 
sive stars, in general, one expects a typical supernovae to 
be much brighter than the star, suggesting that a modest 
contribution to the IRB could come from supernovae, when 
compared to the fractional contribution from stellar emission 
alone. However, though supernovae are in fact brighter than 



an average star, the duration over which most optical fiux is 
emitted is relatively smaller. We show that this results in a 
small but non-negligible contribution to the IRB from super- 
novae associated with massive Pop HI stars. While the back- 
ground contribution is smaller, these supernovae are likely 
to be detectable directly in deep IR images above 1 /xm down 
to AB-magnitude limits at the level of 26 and fainter. 

In the next section, we discuss the contribution to IRB 
from the first generation of stars, supernovae, and mini- 
quasars under the assumption that the early stellar pop- 
ulation is dominated by very massive stars. In the case of 
supernovae, we will suggest that, while the relative contri- 
bution to the IRB may be lower, individual detections may 
be possible with deep-IR imaging. We also suggest that the 
IR background itself may be used to constrain the presence 
of early mini-quasars. Throughout the paper, we make use 
of a ACDM cosmological model consistent with current data 
(e.g., Tcgmark ct al. 2004) with Qm = l-f^A = 0.3, h = 0.7, 
Qi = 0.04, the spectral index of the primordial power n = 1, 
and a normalization for the matter power spectrum with 
0-8 = 0.9. 



2 FIRST-SOURCES IN THE IR 
BACKGROUND 

Our calculations of the first-source contribution to the IRB 
follow previous calculations in the literature related to the 
stellar contribution to the IRB (Santos et al. 2002; Sal- 
vatcrra & Fcrrara 2003; Cooray et al. 2004; Kashlinsky et al. 
2004). These studies involve two basic ingredients: (1) the 
rate at which the volume density of first-sources evolves and 
(2) the flux spectrum. We extend these calculations to su- 
pernovae and miniquasars, which have been so far ignored 
in previous estimates, under the assumption that (1) the 
supernova rate in the early universe traces that of the star 
formation and (2) the formation rate of miniquasars is sim- 
ply related to the formation rate of seed black holes, taken 
to be proportional to the star-formation, and to a model de- 
scription for accretion. We use a simple model to describe 
the spectrum and the light curve of supernovae following 
results from Heger et al. (2002). 

For miniquasars, wo make use of the average spectrum 
of quasars derived by Sazonov ct al. (2004) using fits to ob- 
servational data over the range of wavelengths correspond- 
ing to X-rays and optical and, thereafter, using models to 
describe the emission, such as in the IR band. While the 
spectral shape is well defined in the average sense for ob- 
served luminous quasars at redshifts less than ~ 6, it is 
unclear to which extent this spectrum cam be applied to 
miniquasars at redshifts greater than 10 and with masses 
in the range of lO'^ Mq to lO"* Mq. At energies below 13.6 
cV, the Sazonov ct al. (2004) spectrum, lyFi,, scales as u^'^ '' 
though expectations for miniquasars are that the spectrum 
is harder such that vF^, is a constant (Madau et al. 2004; Di- 
jkstra, Haiman & Loeb 2004). Regardless of the exact shape 
of the spectrum, there is also an uncertainty with respect 
to the overall normalization, or the total luminosity of a 
given miniquasar. Here, we make the assumption that mini- 
quasars radiate efficiently at the maximum allowed by the 
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Figure 1. Top panel: The global star formation rate density at 
redshifts greater than 10 calculated based on the Press-Schechter 
(PS) mass function with a minimum virial temperature of 400 
K (top line) and 10* K (bottom line). Bottom panel: The baryon 
fraction in black holes as a function of redshift, under the assump- 
tion that the star-formation history related to the seed population 
follows the two curves shown in the top panel. We have assumed 
a 5% of the mass in stars are converted to seed black holes. The 
three curves in each of the star-formation models follow a model 
description related to the duty cycle of accretion growth following 
Ricotti & Ostriker (2003) with parameters shown on the figure 
(see Section 2.1 for a discussion). 



Eddington- limit. In making this choice, we are guided by 
observations of quasars out to redshifts of 6 which indicate 
maximal emission at the Eddington- limit. These quasars are 
more massive and more luminous than the miniquasar coun- 
terparts at high redshifts. It is likely that the miniquasar 
emission is submaximal. The limits we derive here on the 
fraction on miniquasars present during the reionization era 
should, thus, only be taken as a conservative limit on the 
low end given the current state of knowledge on the IR back- 
ground. 



Figure 2. The flux spectrum, uF^,, of a 100 Mq star at a 
redshift of 15 as observed today as a function of the observed 
wavelength. In addition to the stellar spectrum (solid line), we 
also show the nebular Lyman-Q emission (long-dashed line) and 
free-free (dashed line) of the ionized Hll region surrounding the 
star, when the UV end of the stellar-spectrum (thin solid line) 
is used to ionized the surrounding medium. The spectrum fol- 
lows from model calculations by Santos et al. (2002). We assume 
all ionizing photons are absorbed by the nebula. For comparison, 
we also show the spectrum of a miniquasar of the same mass 
{Mbh sood = IOOMq) at a redshift of 15, as viewed today (solid 
line labeled 'miniquasar'). The spectrum is computed using the 
template derived in Sazonov et al. (2004). The bump around ^ 
lOO/xm shows a possible contribution from thermal reradiation of 
UV photons by dust. While the overall flux from a miniquasar is 
similar to that of a star, we argue that the contribution to IRB 
from miniquasars could dominate stars because of their larger 
number density; quasars live longer than massive stars. 



2.1 Formation rates and background flux 

In order to calculate the formation rate of first stars and 
the subsequent supernovae or seed black holes, we use an 
analytical description for halo formation following the Press- 
Schechter formalism (Press & Schechter 1974) . We write the 
star-formation rate as 

, , , rib d f°° , , , „ dn , 

^^(^^^^a^d^i^'^^^dM' 

where Mmin is the minimum mass of halos in which gas can 
cool, and we take the star-formation efficiency rj = 0.4. We 
consider two cases for Mmin with a minimum virial temper- 
ature of 400 K, involving molecular hydrogen (e.g., Tegmark 
et al. 1997), and lO^K involving atomic hydrogen line cooling 
(Barkana & Loeb 2001). The respective star-formation rates 
are shown in Fig. 1 (top panel). We mention that, although 
more detailed modeling of early star formation is available 
(e.g. Hernquist & Springel 2003; Yoshida et al. 2003), the 
overall feedback effects from mini-quasars on star-formation 
are still uncertain (see, e.g. Machacek, Bryan, Abel 2003). 
We thus prefer using the model described above for the sake 
of simplicity. Note that the scenario with a minimum tem- 



© 0000 RAS, MNRAS 000, 000-000 



4 Cooray & Yoshida 



perature T^h = 400K likely provides a possible maximal 
star-formation rate. 

Given the rate of formation of stars, we can write the 
emissivity per comoving unit volume at a certain wavelength 
A as a function of redshift z as 

jUz) = ^ {U,^} M^) , (2) 

where 1^ is the luminosity per source mass as a function 
of frequency and fago is the lifetime over which this flux is 
emitted; for a single source at a redshift of z with a source 
mass Ms, such as a first-star, the total Luminosity, as a 
function of frequency, is Lv — IvMs and the observed flux 
today is F„ = Lv{l + z)/[4Tvdj^], where di, is the cosmological 
luminosity distance out to a redshift of z. 

The cumulative background is obtained by integrating 
the emissivity over redshift, yielding the specific intensity Iv 

VobJ. = c[ dz^u{z)^rr^, (3) 
Jo d2; 1 + z 

where Uobs is the observed frequency, ^{z) = (1 -I- 2:)'i;obs 
is the redshift scaling of the frequency, and r is the proper 
distance such that dr/dz = + z)H{z)] when the expan- 
sion rate for adiabatic cold dark matter cosmological models 
with a cosmological constant and a flat space-time geome- 
try is H'^{z) = Ho[flm{l + z)^ + Qa]- In terms of the in- 
dividual source fluxes, and their comoving number density 
na{z), the specific intensity of the background can also be 
written as Vohslv — J dz (dV/d^ldz) h'{z)F,^{z)ns{z), where 
ns{z) — (tagc) Tp*{z)/Ms and {dV/dildz) is the cosmological 
volume element given by d\/{\ + z)^ dr/dz, in terms of the 
luminosity distance. This expression and equation (|3J are 
equivalent. 

We use the same formalism to compute the background 
radiation intensities from stars, supernovae and miniquasars. 
We assume that the supernovae rate follows that of the star- 
formation. For miniquasars, we follow Ricotti & Ostriker 
(2003) and write the growth rate of black holes as 

t^WA^z) = (o^c{z) -I- uj^cca{z) , (4) 

where the two terms represent the growth of black holes 
via accretion, Wac, and the formation rate of seed black 
holes, ajsccd. The over-dot represents the derivative with re- 
spect to proper time. We ignore black hole removal pro- 
cess, such as through ejection from proto-galaxies, and take 
a;seed(z) = /bh^* («) where "(/".(z), the star-formation rate, 
is given by equation (1), and /bh represents the overall frac- 
tion of stars converted to seed black holes. We model growth 
related to accretion as il;ac(z) = /duty(z)wBH(2)/tEdd, where 
Eddington time is set to be a constant at 10* years and 
use the same parametrization in Ricotti & Ostriker (2003) 
for /duty such that /duty(z) = [(1 -I- 2;)/2bh]"* with the 
condition that /duty(z) < 1 based on the definition that 
/duty = ion /(ion + ^off) where ton and toff are time inter- 
vals over which the blackhole accretes and does not accrete, 
respectively. At the low end, following Ricotti & Ostriker 
(2003), we set /duty(-2) > 10~^ to be consistent with both 
the AGN fraction at z ~ 3 and today (Steidel et al. 2002). 
In Fig. 1 (bottom panel), we show ujmi^z) normalized to the 
baryon density and assuming /bh = 0.05 for a variety of 



models related to /duty with parameters as labeled on the 
figure. 

2.2 Source Spectra 

The spectrum of stellar emission is computed as in Santos 
et al. (2002), including the nebular and free-free emission. 
For the miniquasar spectrum, we make use of the model 
of Sazonov et al. (2004). We assume that the Sazanov et 
al.'s template for the average quasars also applies to mini- 
quasars. In Fig. 2, we compare the spectrum of a typical 
star with IOOMq and that of a miniquasar at a redshift of 
15, as viewed today. Here we plot i^obs^l' (in units of W 
m~^ sr~^). The miniquasars are assumed to radiate at the 
Eddington-limit such that the total integrated luminosity 
over the range from X-rays to far-infrared is ~ 1.3 x 10"^* 
ergs s~^ (Mbh/Mq), though, as discussed in Section 2, it is 
likely that the miniquasar emission is submaximal. Also, at 
observed wavelengths below ~ 2 /im, corresponding to red- 
shifted Lyman-Q emission, the miniquasars flux spectrum 
vFv scales as z/""'^. For miniquasars, it could be that the 
spectrum is harder (in terms of energy), such that the ex- 
pectation is that vFv is constant at energies above 13.6 eV 
(Madau et al. 2004; Dijkstra, Haiman & Loeb 2004). This 
will result in a flatter spectrum, as observed today, from op- 
tical/UV to IR wavelengths. As shown in Fig. 2, we assume 
that stellar emission is responsible for reionization of the uni- 
verse; the Lyman-Q emission is related to ionizing photons 
that are first absorbed by the neutral medium and are re- 
emitted during recombinations. In the case of miniquasars, 
we do not consider UV absorption as we have implicitly as- 
sumed that the surrounding medium is already ionized by 
stars preceding the formation of a miniquasar in or near that 
location. We will comment on this later in the discussion. As 
shown in Fig. 2, the level of the incident flux from a typical 
star and a miniquasar is similar. 

The spectrum and light curve for Population III super- 
novae were calculated by Heger et al. (2002). In these cal- 
culations, the relevant IR emission, as observed today, char- 
acterizing the peak of the light curve, is delayed a month or 
so from the initial explosion and the resulting shock break 
out. The latter includes emission at short wavelengths, be- 
low the Lyman limit, in the supernova frame, which is likely 
be absorbed by the intergalactic medium at epochs prior to 
complete reionization. The peak of the light curve extends 
over a month at most. In order to make a reasonable esti- 
mate of the supernovae contribution to the IRB, we use the 
peak emission and its duration only and ignore detailed as- 
pects of the light curve when fluxes drop below more than a 
factor of 10 from the peak emission. Including the full light 
curve only leads to minor corrections at a a few percent 
level. While the peak emission is flat such that \F\ is a con- 
stant in the estimate by Heger et al. (2002), to study any 
departures, we also allow the spectrum to vary as \F\ oc A° 
with a between -0.5 and 0.5, with a = as the fiducial case. 
When varying the spectral shape, we renormalize such that 
the total flux, or the luminosity, remains constant in the rest 
wavelengths between 10^ A and 10* A corresponding to the 
UV and optical regimes. Here, again, we assume that the 
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surrounding medium is reionized by stellar emission before 
the supernova explosion and that the rest UV emission is 
not absorbed. The reionization by supernovae may not be 
important as the time scale over which the rest UV emis- 
sion is expected is significantly smaller than that of a star 
(a month vs. a few million years, respectively). Note that 
the supernova spectrum of a Pop III star is highly uncer- 
tain in terms of spectral shape at the UV end and more 
numerical models are needed to further address the extent 
to which supernovae may be an important source of ioniz- 
ing photons. A more likely scenario involving Pop III super- 
novae is the one described in Oh, Cooray & Kamionkowski 
(2003), where the blast wave propagates to the dense ionized 
IGM and heats the electrons to a substantial temperature 
which subsequently cools via inverse Compton-scattering off 
of the cosmic microwave background, similar to the same ef- 
fect related to hot electrons in galaxy clusters (the Sunyaev- 
Zel'dovich effect; Sunyaev & Zel'dovich 1980). 

2.3 Infrared background 

We summarize our results in Fig. 3, where we show the 
contributions to total IRB light from stars, supernovae and 
miniquasars. In the case of stars, we divide the contribution 
to the three main components; direct emission from stars, 
Lyman-a recombination radiation from the ionized patch, 
and the free-free component related to electron-ion scatter- 
ings within each ionized patch surrounding these stars. Our 
calculations related to the IRB light due to stars arc con- 
sistent with those in Santos et al. (2002) and Salvaterra & 
Ferrara (2003), except that here we have assumed a slightly 
different redshift limit to the first generation of stars at a 
value of 10 instead of values such as 7 and 8 used in those 
calculations that are aimed at explaining the total missing 
content. 

Our results suggest that, while the emission associated 
with stars can explain most of the IRB intensities, the frac- 
tional contribution from supernovae is an order of magnitude 
smaller than the stellar emission contribution. The relative 
difference between stars and supernovae is easily understood 
based on the typical flux from a star vs. a supernovae and 
the ages over which these fluxes are emitted. For example, 
typical ~ 200 Mq star at a redshift of 20 has a flux, ig- 
noring Lyman-a emission, of order ~ 10""^^ ergs cm~^ s~^ 
Hz~^ (~ 10~^ njy or tuab ~ 39), while the peak flux of 
the supernovae at the same redshift is at the level of ^ few 
times 10^"^" ergs cm^^ s^^ Hz^^ {^^ 1 /ijy or rriAB ^ 26). 
Even if all the stars at redshifts greater than 10 die as super- 
novae, though the instant flux is greater, the total output of 
{Fvta.ge) from supernovae is smaller than stars as the stellar 
emission lasts over a million years while the peak flux from 
supernovae lasts just over a month or so. Thus, the apparent 
high ratio of the peak luminosity is largely compensated by 
the age difference so that the typical ratio of a supernova to 
stellar flux remains at most at the level of ~ 0.1. 

Note that the stellar contribution to the IRB also in- 
cludes a substantial fraction from the Lyman-a emission of 
absorbed ionizing UV photons during subsequent recombi- 
nations. We have not included such a contribution from in- 



dividual Pop III supernovae because of uncertainty in the 
supernovae spectrum at the UV end. Note also that re- 
gions surrounding supernovae is likely to have been ionized 
by the progenitor stars. Even if the IGM between proto- 
halos are partly ionized, the resulting Lyman-a emission 
is expected to be suppressed relative to the case of stel- 
lar sources. This can be understood from the fact that the 
Lyman-a spectrum, as a function of frequency, can be writ- 
ten as lhya{i^) — qh h'-Lya4>{^)^ where 4>{v) is the fit to the 
Lyman-a profile of Locb & Rybicki (1999) by Santos ct al. 
(2002) and the production rate of Lyman-a photons from 
recombinations is g ~ 0.6(tage)as(T')ne-R, where Tie is the 
electron density of the surrounding medium, aniT) is the 
case-B recombination coefficient for HI as a function of the 
electron temperature T, and R is the ionizing photon pro- 
duction rate (see Santos et al. 2002 for details). To estimate 
R, we assume that the supernovae spectrum is flat (in rest 
vF,j below the Lyman-limit) and given the uncertainty in 
spectral shape, we take the production rate to be that asso- 
ciated with ionizing photons at the wavelength correspond- 
ing to the Lyman-limit. The model related to stars is shown 
in Fig. 2, where the unabsorbed spectrum peaks at a wave- 
length of ~ 300 A (rest frame) . The difference in flux, as sug- 
gested above, translates to a stcllar-to-supernovae ionizing 
photon production rate difference of ~ 10^'' to 10~®. This 
should then be compared to the ratio of ages, which is again, 
~ lO'^ between stars and supernovEie. Furthermore, in the 
case of supernovae, the Lyman-a radiation is more likely to 
be associated with tenuous IGM with density ~ 10~^(l-|-2;)^ 
cm~^ compared to dense nebulae surrounding stars (with 
density 10"* cm~^). This lowers the Lyman-a flux by another 
factor of 10® relative to that of stars, resulting in an overall 
reduction of lO^'"* to 10^^ in the Lyman-a contribution from 
stars to supernovae. Thus, it is unlikely that the difference 
between predicted supernovae background and the missing 
IRB flux can be reconciled with the expected increase as- 
sociated with Lyman-a emission if the UV emission from 
supernovae are absorbed during the reionization process. 

Intriguingly, redshifted light from early miniquasars can 
explain the IR excess at wavelengths around 1 jj-m, as shown 
in Fig 3(c). In fact, if all stars collapse to become seed black 
holes, we find that the fiux from miniquasars that grow down 
to a ~ 10 will exceed the measured IR fluxes by an order of 
magnitude or more. A fractional contribution of ~ 5% to 
10% of initial stellar mass to seed black holes appears con- 
sistent with observations. This is based on our assumption 
that the miniquasars are radiating efficiently at the maxi- 
mum given by the Eddington-limit. As discussed at the be- 
ginning of Section 2, the extent to which such an assump- 
tion, based on observations of luminous quasars at redshifts 
less than 6, applies to high-redshift miniquasars is unclear. 
If the radiation from miniquasars, with masses between 10^ 
Mq and 10* Mq, is submaximal, the derived limit on the 
fractional contribution from ~ 5% to 10% can be increased. 
Thus our result should be considered as a weak constraint. 

Since the typical stellar age is ~ 2 x 10® years (Bromm, 
Kudritzki & Loeb 2001), at a redshift of around 10, using the 
highest star- formation rate in Fig. 1 (top panel), one finds a 
typical density of order a few times 10® Mq Mpc~^ in stars. 
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Figure 3. The cosmic infrared background due to Population III stars (a), first supernovae (b), and miniquasars (c). The data points 
and upper limits are related to unexplained IRB intensities from measurements and calculations in the literature. For reference, we show 
the same observational data as the ones plotted in Fig. 6 of Santos et al. (2002; the measurements at 1.25 fim and 2.2 fim from Cambresy 
et al. 2001; upper limit at 1.25 fim, a second measurement at 2.2 /im, and a upper limit at 3.5 fj-m from Wright & Johnson (2001); upper 
limits at 0.6 /im and 0.8 fira from Bernstein, Freedman &; Madore 2002). Note that these background measurements are corrected for the 
known contribution from Galactic stars and galaxies down to magnitude limits of order 25 (see, Santos et al. 2002 for further details). 
In the case of stellar contribution to the IRB shown in Panel (a), we show the contribution from stellar emission alone (dotted line), 
associated Lyman-a emission due to recombinations in the ionized patches surrounding individual stars (long-dashed line), the free- free 
emission due to electron/ion scattering within individual ionized patches (dashed line), and the sum of these three contributions (thick 
solid line). These curves make use of the star formation rate (SFR) density based on the PS mass function with a minimum temperature 
for collapse of 400 K; for illustration, the bottom dot-dashed curve is the total contribution to the IR background in the case where 
minimum temperature is set at 10* K. In the case of supernovae, in Panel (b), the three curves assume flux spectra for supernovae 
emission, XF\ oc A° with a = (solid line) ,-0.5 (dashed line), 0.5 (long-dashed line); The curves are normalized such that the total flux 
or total luminosity is same in the UV to optical bands (see Section 2.2 for details). The top curves assume PS mass function to calculate 
the SFR with Tvir=400K, but the case with Tvir=10''K, for a = 0, is shown as a dot-dashed line. In the case of miniquasars in Panel 
(c), we consider three models for the accretion duty cycle as plotted in Fig.l (bottom panel) with parameters given in that figure and 
plotted as top three curves with thick lines when Tvir=400 K. The middle dot-dashed curve is for the case corresponding to the solid line 
but with Tvir=10''K. These curves assume a seed black hole mass fraction, /bhi of 5%. For reference, we also show the range implied 
by the solid line and the dot-dashed line, with thin bottom lines, when this fraction is reduced to 1%. Note that we have assumed the 
case that no UV photons are absorbed by the IGM surrounding these miniquasars. The solid curve labeled '1 a detection limit' is the 
extent to which the optical to near IR background can be established, over the wavelength regime indicated by the curve, using a low 
resolution spectrometer on a proposed rocket experiment. 



From the bottom panel of Fig. 1, one finds that the density 
of miniquasars is typically ~ 10* Mq Mpc""^. This apparent 
increase in miniquasar density, however, is compensated by 
the fact that the miniquasars flux, at relevant wavelengths, 
is roughly an order of magnitude smaller than that of a 
typical very massive star (see Fig. 2). The exact numerical 
calculation reveals that the miniquasars do dominate the 
background if the fraction of stars converted to seed black 
holes is large. In other words, the IRB offers a useful way 
to constrain the presence and abundance of miniquasars at 
z > 10; under our assumption, an initial mass fraction of 
black holes, relative to stars, of order 10% and below seems 
consistent with current observations. 

This conclusion may indeed be strengthened if ionizing 
photons from miniquasars are absorbed by the neutral IGM; 
if absorbed, we expect some fraction of the energy to be rera- 
diated in the form of the Lyman-a line and to increase the 
IRB at the corresponding frequencies such that the upper 
bound on the black hole mass fration is further reduced. 
On the other hand, as discussed in the case for supernovae 
above, we do not expect the Lyman-a radiation to be signif- 
icant as long as the dense nebula surrounding quasar is fully 
ionized already by the stellar emission prior to the forma- 
tion of the miniquasar. The miniquasar emission can ionize 



the IGM beyond the nebula, but due to the lower density, 
recombinations are unlikely to be signiflcant such that the 
Lyman-a emission is suppressed relative to the case where 
stars form and ionize dense nebulae surrounding such stars 
initially. 

Our suggestion that the seed mass fraction must be 
below 10%, assuming efficient radiation at the Eddington- 
limit, also suggests that the universe cannot be reionized 
by early quasars alone. For reference, in Fig. 3, we show 
the model in which the UV photons from miniquasars are 
absorbed by neutral IGM during the reionization process 
(thick dot-dashed line in the figure) . Integrating the UV light 
below the hydrogen Lyman-limit, we find a ratio of num- 
ber density of photons relative to mean density of baryons 
at the level of 0.3 at z ^ 10; to fully reionize, one would 
expect this ratio to be at least unity or as high as 25 
as the recombinations are important at high redshifts. In 
Oh, Cooray & Kamionkowski (2003), it was estimated that 
25± 12(/osc/0.3)"^(re/0.12)(Cii/4) photons per baryon are 
needed to explain an optical depth to reionization of 0.12 
when the clumping factor of ionized gas Cu is 4 and the es- 
cape fraction of ionized photons is 0.3, consistent with other 
estimates (e.g., Dijkstra, Haiman & Loeb 2004). Note also 
that ionization by a combination of hard X-rays and sec- 
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ondary electrons does not lead to complete ionization, as 
argued by Madau et al. (2004). 

Our constraint on the total mass fraction of miniquasar- 
seed blackholes is independent of similar suggestions in the 
literature, for example, based on the unresolved X-ray back- 
ground (Dijkstra, Haiman & Loeb 2004). These authors ar- 
gue that the X-ray background would be violated by a pop- 
ulation of miniquasars if they are responsible for complete 
reionization. If we consider a scenario where miniquasars ex- 
ist with stars, with a seed miniquasar mass fraction of 10%, 
we find that the X-ray background constraint at 1 keV is not 
violated by this miniquasar population even in the extreme 
case where the spectrum is normalized to the Eddington- 
limit. Interestingly, the level of the IRB background from 
the same population is consistent with current observations. 
Given the uncertainties in the available observations and in 
our analytical model, we are not yet able to exactly establish 
the allowed seed mass fraction. 

Although current observations of the IRB are still un- 
certain, one can improve the limit related to the mass frac- 
tion of seed black holes further by future background mea^ 
surements. In Fig. 3(c), we show the Icr noise, or detec- 
tion, level related to a plarmed low-resolution spectroscopic 
background measurement on a rocket experiment (T. Mat- 
sumoto, private communication) for IRB fluctuation studies 
(e.g., Magliocchetti et al. 2003; Cooray et al. 2004). The de- 
tection limit comes from the addition of 400 spectra over 
a 30 second integration time at an altitude above 300 km 
(which removes the air-glow) and after a foreground source 
subtraction down to 14, 15.9 and 16.3 magnitudes in I-, J- 
and K-bands. The background can be measured to a level of 
1 nW m~^ sr~^, whereas the limiting factor in these mea- 
surements will be related to the removal of the zodiacal light 
(e.g., Cambresy et al. 2001). In this case of the rocket exper- 
iment, this removal is achieved with data taken as a func- 
tion of the latitude such that, when combined with multifre- 
quency aspect, removal down to the instrumental noise level 
is expected. The proposed measurement would help estab- 
lish the excess above Galactic stars and foreground sources 
better than what is currently known and, in return, limit the 
extent to which high-z sources, either in the form of stars 
and/or miniquasars, may be present. 

2.4 Far-Infrared background 

One may consider the possibility whether the constraint 
from the IRB can be relaxed if the rest UV and optical emis- 
sion from miniquasars is heavily absorbed by dust. First, 
we note that the extinction by inter-galactic dust is unim- 
portant at ^; > 10 even if nearly all the first stars die as 
pair- instability supernovae (see Yoshida et al. 2004). If the 
inter-stcUar absorption due to dust within (proto-)galaxies 
that are hosting miniquasars is significant, the reradiation by 
thermal dust at far-IR wavelengths at ~ 100/im and above 
can still be used to put limits on their abundance. Assum- 
ing the same spectrum as shown in Fig. 2 for quasars, which 
allows for significant dust absorption, we predict a 850 /itm 
background of ~ 50 nW m~^ sr~^, for /bh = 5% while the 
observed background is below a few nW m~^ sr~^ (Hauser 



& Dwek 2001). This suggests that even for the 5% seed mass 
case, which may explain the missing IRB, absorption by dust 
cannot be as significant as for low redshift luminous quasars. 
Keeping the same seed mass fraction, one can refine the ex- 
tent to which dust absorption is significant by studying the 
unresolved background at these far-IR wavelengths. While 
~ 40% of the 850 /^m and 350 /^m backgrounds have now 
been resolved, the counts axe such that a small extrapolar 
tion to a slightly lower flux density resolves the whole back- 
ground (Borys et al. 2003). Unlike the near-IR background 
where surface brightness converges down to AB magnitudes 
of ~ 25, as observed counts go deeper, such a convergence is 
not found at far-IR wavelengths. Thus, it'll be useful to re- 
turn to this aspect in the future when source counts are well 
deflned. As an early conclusion, we are forced to consider the 
possibility that the far-IR background can be fully explained 
with resolved source counts alone. This, when combined with 
IRB measurements, suggests that a substantial population 
of "dusty" miniquasars is unlikely to be present at z > 10. 

3 SUMMARY 

The cosmic infrared background (IRB) at wavelengths be- 
tween 1 /im and 3 /im provides a useful probe of the global 
star-formation rate prior to and during the reionization. We 
have studied the relative contributions from first stars, su- 
pernovae and early miniquasars to the infrared background 
(IRB). In addition to massive stars, the IRB flux may be 
dominated by a population of miniquasars, especially if the 
first generation of stars are very massive and the end prod- 
uct of the stars lead to a substantial population of black 
holes that will effectively radiate as miniquasars by matter 
accretion. 

We use the Press-Sclicchtcr formalism that describes 
dark halo formation, and combine it with a simple star- 
formation model to calculate the cosmic star-formation his- 
tory, the supernova rate, and the blackliolo formation rate at 
high redshift. We follow previous calculations of the source 
spectrum for stars, while for Population III supernovae, we 
use a simple description of the peak emission and its du- 
ration. We describe the flux spectrum of miniquasars using 
an average spectrum as observed for luminous quasars at 
redshifts less than 5 based on observed data and models 
by Sazonov et al. (2004). We fix the freedom related to the 
overall normalization of this spectrum by assuming that the 
miniquasars radiate efficiently at the maximum allowed by 
the Eddington- limit. Under this assumption, we find that a 
mass fraction of seed black holes, relative to Pop III stars, at 
the level of ~ 10% and below is consistent with observations. 
This upper limit on the seed mass fraction can naturally be 
extended to a higher value if miniquasars radiate below the 
Eddington- limit, so the limit we derive should be considered 
as on the lower side rather than a strict upper limit. Future 
spectroscopic background measurements exploiting rockets 
will improve the observational uncertainties and hence the 
constraints on the high-redshift miniquasar population. 

While the integrated light from Population III super- 
novae to the IRB is subdominant, one can potentially con- 
sider the possibility of directly detecting individual super- 
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novae in high resolution deep IR imaging data. Given the 
star formation rate density and the expected flux, we es- 
timate the surface density of supernovae to be order few 
tens per sqr. degree over a year, if the star-formation rate 
is close to 0.5 M© yr~^ Mpc~^ at z ~ 10, and assuming 
that all stars are massive (> 100 Mq). For the high redshift 
supernovae, the typical AB magnitudes are at the level of 
26 in 1.5 /im (Heger et al. 2002). These are within reach 
of deep IR imaging with existing large telescopes and from 
space and, in the long term, clearly within reach of missions 
such as the James Webb Space Telescope (JWST), which 
can detect point sources down to a magnitude limit of 31 
around the same wavelengths in a 10* sec exposure in its 16 
arcmin^ field of view ^. 
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